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Confining a binary mixture, one can profoundly alter its miscibility behavior. The 
qualitative features of miscibility in confined geometry are rather universal and shared 
by polymer mixtures as well as small molecules, but the unmixing transition in the bulk 
and the wetting transition are typically well separated in polymer blends. We study 
the interplay between wetting and miscibility of a symmetric polymer mixture via large 
scale Monte Carlo simulations in the framework of the bond fluctuation model and via 
numerical self-consistent field calculations. The film surfaces interact with the monomers 
via short ranged potentials and the wetting transition of the semi-infinite system is of 
first order. It can be accurately located in the simulations by measuring the surface and 
interface tensions and using Young's equation. 

If both surfaces in a film attract the same component, capillary condensation occurs and 
the critical point is close to the critical point of the bulk. If surfaces attract different com- 
ponents, an interface localization/delocalization occurs which gives rise to phase diagrams 
with two critical points in the vicinity of the pre-wetting critical point of the semi-infinite 
system. The crossover between these two types of phase diagrams as a function of the 
surface field asymmetry is studied. 

We investigate the dependence of the phase diagram on the film thickness A for an- 
tisymmetric surface fields. Upon decreasing the film thickness the two critical points 
approach the symmetry axis of the phase diagram, and below a certain thickness A tri , 
there remains only a single critical point at symmetric composition. This corresponds 
to a second order interface localisation/delocalisation transition even though the wetting 
transition is of first order. At a specific film thickness, A tr i, tricritical behavior is found. 

The behavior of antisymmetric films is compared with the phase behavior in an anti- 
symmetric double wedge. While the former is the analogon of the wetting transition of 
a planar surface, the latter is the analogon of the filling behavior of a single wedge. We 
present evidence for a second order interface localization/delocalization transition in an 
antisymmetric double wedge and relate its unconventional critical behavior to the pre- 
dictions of Parry et al. (Phys. Rev. Lett 83 (1999) 5535) for wedge filling. The critical 
behavior differs from the Ising universality class and is characterized by strong anisotropic 
fluctuations. 
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1. Introduction 

Confining a binary mixture, one can profoundly alter its miscibility behavior[m2HHllI]. 
The phase behavior of AB mixtures in pores, slits and films has attracted abiding interest 
from both theorists and experimentalists [ E] alike. We study the interplay between 
(pre)wetting and phase behavior by self-consistent field (SCF) theory[ |9j and Monte 
Carlo simulations [ ITU| ITT] IT2] . Particularly, we focus on situations where surfaces attract 
different components of the mixture. 

The qualitative features of the miscibility in confined geometry are rather universal 
and shared by polymer mixtures as well as small molecules. Symmetric binary polymer 
blends are, however, particularly well suited to study the interplay between wetting and 
miscibility: (i) the wetting transition temperature typically is much lower than the critical 
temperature, where demixing occurs in the bulk[^j, (ii) Fluctuations can be controlled 
by the degree of interdigitation[ EH]: the more extended the molecule is, the larger 
is the number of neighbors it interacts with, and the smaller the effect of fluctuations. 
Therefore SCF calculations provide an accurate description for many properties except 
for the ultimate vicinity of critical points. The spatial extension of the molecules also sets 
the length scale of enrichment layers and facilitates experimental investigations. Indeed, 
wetting transitions have been studied in recent experiments^ [7]. (iii) The vapor pressure 
of polymer films is vanishingly small, hence effects of evaporation can be neglected, (iv) 
Polymers tend not to crystallize easily. Therefore, wetting phenomena might not be 
preempted by crystal phases. Likewise there is no roughening transition of the interface 
as it occurs in Ising-like models. 

Using a coarse-grained polymer [ ITSl HI] model for a AB binary melt we locate the first 
order wetting transition, the phase diagram in a symmetric slit pore (symmetric film) [ITU], 
the phase diagram in a thin film where the substrate favor the A-component of the mixture 
with the same strength as the top surface attracts the S-component (antisymmetric film) [ 
E] • Then we discuss the phase behavior in a quadratic pore where two neighboring 
surfaces favor the ^-component and the other two neighboring surfaces favor the B- 
component (antisymmetric double wedge) [ UJ3 [TE] . We conclude with an outlook. 

2. Model and techniques 

We consider a binary polymer blend. Both species - A and B - contain the same 
number N of monomers and have the same spatial extension R e . They are confined into a 
thin film; the bottom substrate (W) might be a silicon wafer, while the other surface might 
be the interface to the vapor (vacuum, V). Depending on the ratio between the interface 
tension jab between the segregated bulk phases and the surface tension javi Ibv of the 
components and the vapor, the upper surface might be rough. The qualitative behavior 
is illustrated in FigHJ If the AB interface tension is comparable to the liquid/vapor 
tension, it "drags" the film surface towards the substrate so as to reduce the length of the 
AB interface. If the liquid/vapor tension exceeds the AB interface tension by about two 
orders of magnitude, however, the surface is almost flat and the situation is equivalent 
to a binary mixture between two hard walls a distance A apart [EI- I n the following we 
shall restrict ourselves to this limit jab "C 7a v or Jbv- 

In the Monte Carlo simulations we use a computationally efficient, coarse-grained lattice 
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lateral coordinate 



Figure 1. Laterally segregated binary 
film. The shape of interfaces is obtained 
by minimizing the effective Hamiltonian 
H = 'JabLab + IasLas + IbsLbs + 
^avLav +1bvLbv at fixed volume of the 
components. Ly denotes the length of 
the interface between substances i and j, 
and jij the corresponding interface ten- 
sion. Jav ~ 1BV = lAS ~ IBS = 0.57AB, 

Ibs = Iab and 'Jav/iab as indicated in 
the key. From [H7j. 



model. The bond fluctuation model[EIJ retains the universal features of polymers - 
connectivity, excluded volume of segments and a thermal interaction which leads to phase 
separation - but ignores details of chemical structure. Effective monomers prevent the 
corners of a unit cell of a 3D cubic lattice from double occupancy. We use chain length 
N = 32 and R e « 17u. Monomers along a chain are connected via bond vectors of 
length 2, \/5, VQ, 3 or y/lO in units of the lattice spacing u. Different monomers repel 
each other by a square well potential of depth e which comprises the nearest 54 neighbors, 
like monomers attract each other. The strength of the repulsion is proportional to the 
Flory-Huggins parameter \ — 5.3e//csT[|TSj- Surfaces are structureless and impenetrable. 
They act on monomers in the two nearest layers (d wa ii = 2) with strength e wa n. 

In the SCF calculations we model the polymers as Gaussian chains[ |HJ El QH|- The 
repulsion between different species is quantified by the Flory-Huggins parameter x- Short- 
ranged interactions of strength Ai and A 2 attract (repel) the A (B) component in the 
vicinity of the surfaces. The total density profile of the film is imposed. It smoothly 
decays to zero at the surfaces in a boundary region of width 0.15i? e . The blend is assumed 
to be incompressible. This standard Gaussian chain model is solved within mean field 
approximation. 

3. Wetting transition 

To accurately locate the wetting transition and calculate the contact angle of macro- 
scopic v4-drops we use Young's equation[[Tn] 7^_bcos0 = Jwb ~ Iwa- Computationally, 
this techniques [ ITUI l2~0j has distinct advantages for locating first order wetting transitions: 
(i) The interface free energy jab and the difference A7 = ^wb — Iwa can be measured 
accurately in separate simulations thereby avoiding the need for huge simulation cells to 
simulate a thick A-layer at the surface in equilibrium with a 5-rich bulk, (ii) Unlike 
observing the dependence of the thickness of the A-layer on temperature or monomer- 
surface attraction, one directly measures free energies. Therefore, we do accurately locate 
the transition, while the instability of the A-rich layer is located between the transition 
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Figure 2. Interface tension ^ab and 
difference in surface tensions A7 as a 
function of inverse temperature e/ksT 
obtained from simulations. Approxima- 
tions for the interface tension jab = 
b^^/x/Q and A7 = 2$c/ waU e wa n = 
e wa n/4 in the strong segregation limit 
also shown. From [EI]- The inset shows 
the dependence of the contact angle on 
e wa n for the two temperatures investi- 
gated in Sec. 5. 



and the mean-field wetting spinodal. (iii) By virtue of the A ^ B symmetry, the differ- 
ence A7 can also be rewritten as the difference A7 = jwb — 1-wb of surface tensions of 
a wall that attract the A-component and a wall that attracts the 5-component. This free 
energy difference can be measured by thermodynamic integration or expanded ensemble 
methods [ HU] • 

The results for our model are presented in Fig|2] From the crossing of 7ab(c) an d 
A7(e) we locate the wetting transition. The fact that curves intersect under a finite angle 
indicates that the wetting transitions are of first order. As we reduce the monomer-surface 
attraction the wetting transition shifts to higher temperatures fc^T/e and become weaker. 
For all transitions studied, however, the wetting transition is of first order. This is also 
corroborated by SCF calculations [ [§1 E], where we find first order wetting transitions for 
T/T c < 0.98. 

If the wetting transition is of first order, then there will be only a small A-rich layer 
in the non-wet state. By virtue of the structural symmetry of the molecules, they loose 
the same amount of entropy as they pack against the surface. The surface free energy 
difference A7 is mainly enthalpic. If we assume that the wetting transition is strongly first 
order, we can neglect the microscopic enrichment layer at the surface A7 = 2e wa n<i wa n ( 3>, 
where d wa \\ = 2 denotes the range of the monomer- surface interaction and $ = 1/16 the 
monomer number density. Using the expression for the interface tension ^ab = $6a/ x/6 
(b = 3.05: statistical segment length) in the strong segregation limit [ 22 , we obtain: 

Xwct — ^ y bk B T ) 

This is in marked contrast to the value of the Flory-Huggins parameter at the unmixing 
transition in the bulk, Xc = 2/iV ~ 1/T C . As both the interface tension ^ab and the 
difference in surface tension A7 are chain length independent, so is the wetting transition 
temperature. The critical temperature T c of phase separation, however, increases linearly 
with chain length N. Therefore, critical phenomena associated with the bulk unmixing 
and wetting phenomena are well separated. 
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Figure 3. Phase diagram in terms of 
composition and temperature for film 
thickness A « 2.8R e as obtained from 
simulations The arrow marks the wet- 
ting transition temperature. e wa u = 
0.16k B T From [[TT| 



4. Thin films 

4.1. Capillary condensation and interface localization/delocalization 

If the mixture is confined into a film the surface interactions modify the phase behavior. 
As wetting is associated with the growth of an infinitely large enrichment layer, it is 
rounded-off in a thin film[|3]. If the wetting transition is of first order, there will be a 
pre-wetting transition [ [I]: a coexistence between a thin and a thick (but microscopic) 
enrichment layer at a chemical potential which differs from the value at coexistence in the 
bulk. As pre-wetting transitions involve only enrichment layers of finite thickness, they 
might give rise to transitions in thin films. 

First we regard a film with symmetric surfaces[E3 E2] i.e. both surfaces attract the 
A component. The phase diagram as obtained from the simulations is presented in Fig. 
El Compared to the phase behavior in the bulk, the critical point is shifted to lower 
temperatures and larger composition of the species attracted by the surfaces. Moreover, 
the binodal in the vicinity of the critical points exhibit 2D Ising critical behavior in 
contrast to the 3D Ising behavior of the bulk unmixing transition. 

Note the pronounced distortion of the 5-rich binodal in the vicinity of the wetting 
transition. In the S-rich phase there are A-rich layers at the surfaces and the B component 
prevails in the middle of the film. In the vicinity of the wetting transition the thickness 
of the A-enrichment layers grows as we increase the temperature. If we increased the 
film thickness this distortion would evolve into an additional two phase regionf^Ul EH]- 
corresponding to a 5-rich phase with thin and thick A-layers at the surface. This two 
phase region would correspond to the pre-wetting coexistence and it would join the S-rich 
binodal in a triple point. 

The phase diagram of an antisymmetric film is also presented in FigCJl In this case one 
surface attracts the A-component with exactly the same strength than the other surface 
the 5-component. The phase diagram contains two critical points and a triple line[|BllI3 
ITT| . Around the critical temperature of the bulk, enrichment layers gradually form at the 
surfaces and stabilize an AB interface that runs parallel to the surfaces. At the interface 
localization/delocalization transition [|2H 123 EH] this AB interface becomes bound to one 
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Figure 4. (a) Schematic temperature dependence of the effective interface potential in 
a film with antisymmetric surfaces. The temperatures correspond to T < T trip , T trip , 
^trip < T < T c film and T c film . (b) Phase diagram of a mixture, (c) Sketches of typical 
configurations for T > T c film (upper panel), T trip < T < T c film in the miscibility gap 
(middle panel) and T < T trip (lower panel). From [ IT2"] 



of the surfaces. In the case of a first order interface localization/delocalization transition 
this corresponds to a triple point of the phase diagram: an A-rich phase , a I?-rich phase 
and a phase with symmetric composition coexist. 

The behavior can be analyzed qualitatively by looking at the interface potential g(l) 
which describes the interaction between an AB interface and a single surface. If the film is 
thick enough, the interface potential can be constructed as a superposition of the interface 
potentials emerging from each surface. The qualitative behavior in the vicinity of a first 
order wetting transition is depicted in Fig0Ja). Using a double-tangent construction we 
can construct the phase behavior in a thin film. At low temperatures there coexist an 
A-rich phase and a 5-rich phase, in which the AB interface is localized at the surface. 
Upon increasing the temperature, one encounters the triple point. This triple point is the 
thin film analogon of the first order wetting transition. As the film thickness increases 
the triple temperature converges towards the wetting transition temperature of the semi- 
infinite system. Above the triple temperature there are two phase coexistence regions, 
which correspond to thin and thick enrichment layers at the surfaces. This is the analogon 
of the pre-wetting transition in a thin film. 

4.2. The tricritical interface localization/delocalization transition 

If we reduce the film thickness, the interactions emerging from each surface interfere. 
The phenomenological considerations [ explain that this leads to a second order in- 
terface localization/delocalization transition at small film thicknesses. Both regimes are 
separated by a tricritical transition. The scaled distribution functions prove convenient 
to locate the tricritical thickness accurately. To this end we have adjusted the temper- 
ature such that the central peak of the probability distribution of the order parameter 
m ~ <Pa — <Pb is a factor 1.2 higher than the outer peaks. This corresponds to the behavior 
of the universal distribution of the 2D tricritical universality class[|27j. The results for 
various film thicknesses A (in units of the lattice spacing) are presented in FigJSJ For 
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Figure 5. Probability distribution for 
various film thicknesses as indicated in 
the key scaled to unit norm and variance. 
The lateral system size is L = 96. The 
universal distribution of the 2D tricriti- 
cal universality class [|2Zj is presented by 
circles. From [ ITT] 



A < A tr i the valleys between the three peaks are too shallow (cf. Fig|5J), while they are 
too deep for A > A tr j. In the latter case the transition is of first order and our estimate 
tends towards the triple temperature. At A tri « 14 = 0.89R e the distribution of our 
simulations is similar to the universal 2DT distribution, and this has been confirmed for 
larger lateral system sizes [ITT]. 

4.3. Crossover from capillary condensation to interface localization/delocalization 

Realizing strictly (anti) symmetric surface interactions is often difficult in experiments. 
Varying the surface interaction A2iV of the top surface from attracting the A-component to 
attracting the 5-component (while the bottom surface always attracts the A-component 
with fixed strength AiN) in the SCF calculations, we study the crossover from capillary 
condensation for symmetric surfaces to interface localization/delocalization. The depen- 
dence of the phase diagram on the surface interactions within the SCF calculations is 
presented in FigEJ For symmetric surfaces (capillary condensation) the critical point is 
shifted towards lower temperatures [ 4 similar to the simulation result. Of course, the bin- 
odals are parabolic in mean field theory independent from dimensionality. The coexisting 
phases have almost uniform composition across the film and differ in their composition. 
As we reduce the preference of the top surface for species B, the critical point and the 
critical composition tend towards their bulk values (0 = 0.5, 1/xiV = 0.5), i.e. the criti- 
cal temperature increases and the critical composition becomes more symmetric [|Sj. The 
coexistence curve in the 1/xN-Afj, plane approaches the symmetry axis. Upon making 
the top surface attracting the other component B, we gradually change the character of 
the phase transition towards an interface localization/delocalization transition [ EU 125] . 
The critical temperature passes through a maximum and the critical composition through 
a minimum, respectively. For A2A < (surface attracting the 5-component) there are 
enrichment layers of the A-component at the bottom and the .B-component at the top, 
and the two coexisting phases differ in the location of the AB interface which runs par- 
allel to the surfaces. As the preferential interaction of the top surface increases, the 
critical temperature decreases and the critical composition becomes richer in A. When 
the coexistence curve intersects the pre- wetting line of the bottom surface at A/i < 0, 
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Figure 6. (a) Binodals for A = 2.6R e and surface interaction AiN = 0.5 obtained from 
SCF calculations. The surface interaction at the other surface A 2 iV varies as indicated in 
the key. The dashed curve shows the location of the critical points. Filled circles mark 
critical points, open circles/dashed horizontal lines denote three- phase coexistence for 
A2N = —0.3675 and —0.5. (b) Coexistence curves in the xAf-A/i plane. The "quasi-pre- 
wetting" lines for A/x < and A 2 iV = —0.3675 and —0.5 are indistinguishable, because 
they are associated with the pre-wetting behavior of the surface with interaction AiN = 
+0.5. FromflH] 



a triple point forms at which an A-rich phase and two I?-rich phases with a thin and a 
thick v4-enrichment layer coexist. When the bottom surface attracts A with exactly the 
same strength as the top surface B (antisymmetric surfaces), the phase diagram becomes 
symmetric. 

For symmetric surface fields the critical point occurs close to the bulk critical point 
(and converges towards it in the limit of infinite film thickness) while the critical points 
in antisymmetric films are associated with the wetting transition and converge towards 
the pre-wetting critical temperature of the semi-infinite system (if the wetting transition 
is of first order) for 5 — ► 00. In both cases, however, critical points belong to the 2D Ising 
universality class. 

5. Interface localization/delocalization in an antisymmetric double wedge 
5.1. Background 

In the following we consider wetting (or rather filling) in a wedge geometry. Macroscopic 
considerations show that the wedge will be filled with liquid when the contact angle O 
on a planar substrate equals the opening angle a. Intriguingly, Parry and co-workers[ 28j 
predict that the filling of a wedge is related to the strong fluctuation regime of critical 
wetting and that critical filling may even occur if the concomitant wetting transition of 
the planar surface is of first order. Specifically, they predicted the distance l of the 
AB interface from the bottom of a wedge to diverge as Iq ~ (Tf — T)~ f3s with (3 S = 1/4. 
Correlations along the wedge and in the other two directions are characterized by diverging 
correlation lengths £ y ~ (Tf — T)~ Uy and £ x ~ £j_ ~ (Tf — T)~ u± with exponents u y = 3/4 
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Figure 7. Antisymmetric double 
wedge: periodic boundary conditions 
apply along the y-direction and there are 
4 impenetrable surfaces of size L x L y . 
The bottom ones (W\) attract the A- 
component with strength e wa n and the 
top ones (W 2 ) attract the 5-component. 
l denotes the position of the interface 
from one corner. From [El- 



and u± = 1/4, respectively. 

In the following we study a wedge with opening angle a = n/4 of the wedge (c.f. Fig[7J). 
Similar to the study of wetting an antisymmetric geometry is advantageous. Therefore 
we stack two wedges which attract different components on top of each other. This 
antisymmetric double wedge is a pore with quadratic cross-section of size L x L. Let L y 
denote the length of the wedge (c.f. Fig[7J). (i) If we used identical surface fields on all 
four free surfaces the analog of capillary condensation would occur in a wedge, i.e., phase 
coexistence would be shifted away from the bulk coexistence curve and the wetting layers 
would be only metastable (with respect to "wedge condensation"), (ii) As the wetting 
layer grew on all four surfaces in the case of symmetric boundaries, we would need larger 
system sizes to reduce interactions between the wetting layers across the wedge. 

The phase behavior in such an antisymmetric double wedge geometry has been studied 
recently in the framework of an Ising model[E3 E]. When the wetting transition of the 
planar substrate was of first order, the wedge filling was also found to be of first order. 
When the wetting transition was of second order an unconventional scaling behavior was 
observed which is characterized by critical exponents a = 3/4, (3 = 0, and 7 = 5/4. Those 
critical exponents can be related (see below) to the exponents of critical wedge filling, and 
the simulations of the Ising model confirm the predictions of Parry and co-workers [128]. 

In the following we corroborate these findings in the framework of the Ising model by 
our polymer simulations. Moreover, we present evidence for the unconventional second 
order transition in an antisymmetric double wedge even though the wetting transition on 
a planar substrate is of first order. 

We present preliminary simulation data for two temperatures: e/fcgT = 0.025 (T/T c = 
0.58) and e/hgT = 0.05 (T/T c = 0.29). At both temperatures the wetting transitions, 
that occur at appropriate attractive strength e wa n of planar surfaces, are of first order (c.f. 
Fig|2J). In the former case it is weak first order wetting transition, in the latter case it is 
a strong first order transition. 

5.2. First order transition in an antisymmetric double wedge 

At the lower temperature e/ksT = 0.05, the behavior is similar to a first order interface 
localization/delocalization transition. We consider here only the case A/x = where phase 
coexistence in the bulk occurs. This excludes the rather interesting interplay between 




9 




Figure 8. The probability distribu- 
tion of the composition at e/ksT = 0.05 
and system geometry 64 2 x 228 exhibits 
a three peak structure, which is charac- 
teristic of a first order transition. The 
inset shows the dependence of the cu- 
mulant (m 2 )/(\m\) 2 with m ~ <f) — 1/2 
on e wa n for three different system sizes. 



pre- wetting and pre-filling behavior studied in Ref . [ S223 • At large surface interaction 
e wa ii > e ^an' wedge there runs an AB interface along the diagonal which divides the two 
double wedges. This corresponds to the delocalized state. Upon decreasing e wa u (or 
decreasing the temperature) the AB interface becomes localizes in one of the wedges. In 
this case the composition of the double wedge is either A-rich or 5-rich and we define 
as order parameter m = 4>a — 4>b- The two situations are separated by a triple point 
e wa,n' Wedge & t "which the interface can be localized in either of the wedges or delocalized on 
the diagonal. The trimodal probability distribution in the vicinity of the tricritical point is 
presented in Fig|8l In analogy to the case of antisymmetric films we expect this triple point 
in a double wedge to correspond to a first order filling transition. In the inset we show 
the cumulant (m 2 ) / (\m\) 2 . If the transition was of second order, these cumulants would 
depend monotonically on e wa u and would exhibit a common intersection point. This is 
not at all what we observe, and we conclude that the interface localization/delocalization 
transition in the double wedge is of first order at the lower temperature e/fcgT = 0.05. 

5.3. Critical behavior in an antisymmetric double wedge 

Even though the wetting transition on a planar surface at e/ksT = 0.025 is of first order, 
the behavior at the interface localization/delocalization transition in an antisymmetric 
double wedge at high temperature differs from the transition at low temperature. In the 
inset of FigEl (a) we present the dependence of the cumulant on the surface interaction 
strength for various system sizes. The cumulants depend monotonically on e wa n and 
exhibits a common intersection point around e™^ ~ 0.035. In panel (b) we show the 
probability distribution of the composition at this intersection point: the distribution 
is bimodal and the two largest system sizes collapse onto a master curve without any size 
dependent prefactor. Therefore we conclude that the interface localization/delocalization 
transition is of second order. 2 

2 The interface localization/delocalization transition might be of second order in a very thin antisymmetric 
film (c.f. Sec. 4.2) even if the wetting transition is of first order. Therefore still larger system sizes would 
be desirable to confirm this conclusion. We note however, that the thickness of the enrichment layer at 
the first order wetting transition of the planar substrate (e/ksT = 0.0226, eJ^/fceT = 0.04) is only 
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Figure 9. (a) Dependence of the absolute order parameter m = \4>a — 0b I on the surface 
interactions e wa n at e/ksT = 0.05. The insert shows the cumulant. (b) Scaling of the 
probability distribution at e wa n = 0.035 and various system sizes. 



Intriguingly there are also marked differences between this second order transition in an 
antisymmetric double wedge and the second order transition in a thin film which belongs 
to the 2D Ising universality class. In the latter case, only the distribution of the scaled 
order parameter L^^m, where (3 = 1/8 and v = 1 are the critical exponents of the order 
parameter and the correlation length in the 2D Ising universality class, exhibits data 
collapse for different system size. Moreover, we present in FigJH^a) the dependence of the 
absolute value of the magnetization in the vicinity of the transition. Curves for different 
system sizes exhibit a common intersection point which agrees well with the intersection 
point of the cumulants. The analogous curves at a Ising-like transition do not exhibit 
a common intersection point but monotonously converge towards (\m\) ~ \T — T c \@ for 
T < T c and (\m\) = for T >T C upon increasing the system size. 

To relate the critical behavior of the antisymmetric double wedge to the predictions of 
Parry et a/.[|2E], we regard the distance Iq of the AB interface from the corner of one wedge. 
Similar to an antisymmetric film (c.f. Sec. 4.1), we assume that we can approximate the 
distribution in a double wedge by the superposition of the distributions of single wedges 
-Pwedge(^o) via P(l ) ~ P we dg (^o) + -Pwodgc ( V^L - l ) If the two distributions P W ed g e(Zo) and 
-Pwed g e(\/2£ — ^o) do not overlap, the AB interface will be located in either of the two 
wedges and the order parameter will not vanish. If the two distributions overlap, the 
interface fluctuates around the diagonal and the order parameter will be zero. Right at 
the transition the two distributions begin to overlap: 

(^o)+£_l = y/2L— (l ) — (interface localization/delocalization in double wedge) (1) 

where (lo) denotes the mean height in a single wedge and £j_ its fluctuations. Importantly, 
Parry's prediction of (3q = u± in wedges (and also corners [130]) means that the height and 
its fluctuations are of the same order. They diverge as we approach the critical filling 
transition. 



lo ~ 4 <§; 45.2 = L/V2. Therefore we believe that our conclusion is not affected by finite size effects. 
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The height of the interface Iq is related to the order parameter m of the localiza- 
tion/delocalization transition. Therefore we expect the distribution of the order param- 
eter also to be bimodal. As Iq ~ £ ~ L at the transition and the order parameter is 
a function of Iq/L the distribution of the order parameter will exhibit two peaks whose 
positions and widths will not depend on the system size. This is exactly what we observe 
in FigJ^b). Using this observation and the standard finite size scaling assumption at a 
second order phase transitions 

P(m) ~ L^ViL^m, L/£ x , L y /Q ~ L^V{L p ' v ^m, L l/u H, rj) (2) 

where V and V are scaling functions, t = (T — Tf)/Tf denotes the relative distance to the 
filling transition, and r] = L y /L Vv ' Vx = L y /L 3 denotes the generalized aspect ratio, we 
conclude fi/v± = 0. Due to the anisotropy of the fluctuations of the interface along the 
wedge with correlation length £ y and perpendicular to the wedge with correlation length 
the generalized aspect ratio appears as a scaling variable. In our simulations we have 
chosen the system geometry such that r\ remains approximately constant to ensure that 
finite size finite-size rounding in the direction along the wedge and the rounding in the two 
other directions set in simultaneously. 3 Hence, the scaling of the probability distribution 
not only confirms (3 = but also u y = 3v±. 

Knowing the probability distribution of the order parameter we can calculate all its 
moments: 

(m k )=M k (L l ^H,r ] ) (3) 

where M. k are scaling functions. As a special case, we calculate the susceptibility: x — 
L 2 L y (m 2 )/k B T ~ L 2 L y M 2 (L/£ ± ,L y /Q ~ £fo _ t -2u ± - Vy = r7 with 7 = 2l/J _ + ^ = 
5/4. Gratifyingly these values for the exponents comply with the anisotropic hyperscaling 
relationf 31 7 + 2/3 = (d — l)u± + u y . Using thermodynamic scaling 2 — a = 7 + 2/3 we infer 
the critical exponent a = 3/4 for the specific heat. Another consequence of the absence of 
any L-dependent prefactor in Eq.Q is the common intersection of moments of the order 
parameter at the transition. Again this is an agreement with our observation in FigJHJ^a). 
As this intersection involves only the lowest moment of the order parameter it yields 
an accurate estimate of the location of the critical interface localization/delocalization 
transition in an antisymmetric double wedge. 

It is interesting to relate the observation of first and second order interface localiza- 
tion/delocalization transition in a double wedge to the shape of the interface potential. 
Parry et al. predict [|2E] that the filling transition is second order if the interface poten- 
tial between an AB interface and a planar surface does not exhibit a free energy barrier 
between the minimum close to the surface and the behavior at large distances, i.e. if a 
macroscopically thick film is not even metastable. 

In FigElwe present the interface potential obtained from the probability distribution 
of the composition in a simulation of an antisymmetric film at e/ksT = 0.025. In the 
vicinity of the wetting transition the interface potential exhibits a maximum between the 

3 If we kept the ratio L y /L constant 77 — > and the system would exhibit a behavior characteristic of a 
corner. In the limit L fixed but L y the wedge becomes quasi-onedimensional and there is no transition. [ 

on 
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Figure 10. The probability distribution 
of the composition in an antisymmetric 
film with system geometry 48 2 x 96. Val- 
ues of e/ksT and e wa a/k,BT (shown in 
the key) correspond to the wetting tran- 
sition and the filling transition (accord- 
ing to Young's equation). 
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minimum close to the surface and the value at large distances. This fact confirms that the 
wetting transition is of first order. At the smaller value of e wa n, however, there is no such 
maximum within the statistical uncertainty of the Monte Carlo data and, in agreement 
with Parry's predictions, we observe a second order transition in the double wedge. 

6. Summary 

We have investigated the interplay between wetting and phase separation of incompress- 
ible binary mixtures confined in thin films and wedges. In our polymer model the wetting 
transition if of first order and we can accurately locate it via Young's equation. The con- 
comitant pre- wetting behavior modifies the phase boundaries in thin films. If both surfaces 
attract the same component, capillary condensation occurs and the critical point is close 
to the critical unmixing transition in the bulk. If one surfaces attracts the A-component 
but the other attracts the 5-component an interface localization/delocalization transi- 
tion occurs. In this case there are two critical points which correspond to the pre-wetting 
critical points at each surface. If the film thickness is very small, however, the interface 
localization/delocalization transition might be of second order even if the wetting transi- 
tion is of first order. The critical points in a thin film are characterized by Ising critical 
behavior. 

In analogy to the interface localization/delocalization in an antisymmetric film, we 
have studied the transition in an antisymmetric double wedge and we relate the phase 
behavior to the filling transition in a single wedge. Importantly we present evidence 
that the analog of critical filling in an antisymmetric double wedge geometry gives rise 
to unconventional critical behavior characterized by an order parameter exponent (3 = 
and strong anisotropic fluctuations[ H^j. We can relate the critical exponents to the 
predictions of Parry et a/.[|2E| on critical filling. In agreement with those predictions the 
filling transition can be critical even though the wetting transition on a planar substrate 
is of first order. This is practically important because there is no experimental realization 
of critical wetting on a solid substrate. Our findings suggests the polymer blends might 
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be promising candidates to explore the filling behavior experimentally. 
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